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Abstract 

We investigate the density profile in a driven diffusive system caused by a 
plane particle source perpendicular to the driving force. Focussing on the 
case of critical bulk density c we use a field theoretic renormalization group 
approach to calculate the density c{z) as a function of the distance from the 
particle source at first order in e = 2 — d (d: spatial dimension). For d = 1 
we find reasonable agreement with the exact solution recently obtained for 
the asymmetric exclusion model. Logarithmic corrections to the mean field 
profile are computed for d = 2 with the result c{z) — c ^ z" -i(ln(z))2/3 for 
z oo. 
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I. INTRODUCTION 



It is well-known that in thermodynamic systems with longe-range correlations boundaries 
have a considerable influence on physical quantities even at macroscopic distances from the 
surface. During recent decades this effect has been extensively investigated in the context of 
surface critical phenomena (a review is given in Ref. |T|). It has been shown that the critical 
behavior near a boundary is governed by universal scaling laws with new critical exponents 
which cannot be expressed in terms of bulk exponents. The renormalization group has 
proved to be a useful method for the classiflcation of both static and dynamic |^-|^ 
surface universality classes. 

While in equilibrium systems long-range correlations occur if the thermodynamic param- 
eters approach a critical point, they seem to be the rule in non-equilibrium systems with 



conserved dynamics P,pl)|. In this paper we study the diffusion of particles subject to a driv- 
ing force in a system with open boundaries. Here the particle conservation in conjunction 
with the deviation from detailed balance leads to long-range correlations and anomalous 
long-time behavior |Tl| even at temperatures T ^ Tc above the critical point (for a review 
see [|10|)- Especially interesting from a physical viewpoint are surfaces which act as particle 
reservoirs and thus break the conservation law. In the case of boundaries perpendicular to 
the driving force which we consider in the present paper a particle reservoir is necessary to 
maintain a non- vanishing steady current. 

A simple microscopic realization of a driven diffusive system (DDS) is a lattice gas 
with hard core repulsion and nearest neighbor hopping. In a homogeneous state of density 
(particles per lattice site) c the external field (favoring jumps in the positive z direction) 
produces a steady particle current j(c). Due to the excluded volume constraint the current 
vanishes if every site is occupied by a particle, i.e. j(l) = 0, while j(c) ~ c for c ^ 1. In 
a work by Krug on boundary- induced phase transitions the general form of the density 
profile has been discussed. It has been shown that the bulk density c satisfies the following 
maximum current principle: If the boundary at 2; = is in contact with a particle reservoir 
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of density c(0) and every particle that reaches the boundary at z = L leaves the system, i.e. 
c(L) = 0, then the current is maximized in the sense that 

j(c) = max{j(c) I < c < c(0)} (1) 

for L — s> oo. As a direct consequence of this principle the bulk density is at the maximum 
point c* of the function j(c) and the profile decays only algebraically from boundary- to 
bulk- value if c(0) > c'' > c{L). 

Recently the exact density profile in one dimension has been calculated for arbitrary 



boundary conditions < c(0),c(L) < 1 ||l^p!5| . These works confirm and generalize a large 
part of the results obtained in Ref. 

Until now no exact solutions have been found for more complicated problems like driven 
diffusion in higher dimensional systems, at a critical point |jl6|, or in a medium with quenched 



disorder ||T^. In these cases the field theoretic approach is useful to obtain systematic 
approximations for density profiles or correlation and response functions. In the present 
paper we use the exact solution of the one- dimensional asymmetric exclusion model found 



m 



I^JTH] to test the accuracy of the renormalization-group improved perturbation theory. 
We also extend the analysis to two-dimensional DDS. 

In the next section we present the semi-infinite extension of the continuum model for DDS 



introduced in Ref. [0. In Sect. |IT| the boundary conditions are discussed and the density 



profile is calculated in a mean field approximation. Having introduced the renormalization 



we 



factors which are required to obtain a well-defined renormalized field theory in Sect. |V 
calculate the Gaussian fluctuations around the mean field profile. We use the renormalization 
group to compute the universal scaling function for the profile at first order in e = 2 — 
and compare our result with the exact one-dimensional profile. For = 2 we obtain the 
logarithmic corrections to the mean field solution. Sect. |V] contains a discussion of our 
findings and an outlook. In Appendix ^ some technical details of the calculation of the 
surface divergencies at one-loop order are given. In order to compare our perturbative result 
with the exact solution we have to take the continuum limit of the exact profile. This is 



done in Appendix |B[ 



II. THE MODEL 



Some time ago, van Beijeren, Kutner, and Spohn [jT8| introduced a continuum model 
for the diffusion of particles subject to a driving force. The anomalous long time behavior 
of this driven diffusive system (DDS) was studied by Schmittmann and one of us |[Tl]] by 
renormalization group methods. The equation of motion for the particle density c(r, t) in 
this model is given by the continuity equation 

^c(r,t) + Vj(r,t) = 0. (2) 

Here the current j(r, t) consists of a diffusive part, a contribution caused by the driving force 
E, and a part }ii{r,t) which is assumed to summarize the effects of the fast microscopic 
degrees of freedom: 

j(r, t) = -DVc(r, t) + K(c(r, t))E + hir, t). (3) 

This form of the current can — in principle — be derived from the microscopic dynamics by a 
suitable coarse graining. Since the external field E introduces an anisotropy into the system 
the coarse graining in general gives rise to an anisotropic matrix of diffusion constants D. 

Expanding the conductivity /t(c) in the deviation s = c—c of the density from its uniform 
average (bulk-) density c, we have 

k(c) = k(c) + k{c)s + ^k"{c)s'^ + .... (4) 

Neglecting higher order terms in this expansion one obtains the Langevin equation 

dts{r,t) = X{A^+pdl)s{r,t) 

+XdJfs{r,t) + lgsiv,tY)+ar,t), (5) 
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where / oc —Ek,'{c) and g oc —Ek,"{c), and the indices || and _L distinguish spatial directions 
parallel ('longitudinal') and perpendicular ('transverse') to the driving force. The Langevin 



force ( = — Vj_R is assumed to be Gaussian with zero mean and the correlations (after a 
suitable rescaling of s) 

(C(r, t)C(r', t')) = -2A( + adl)6{r - r')6{t - t'). (6) 



It was shown in |TT| that the model defined by Eqs. is complete in the renormalization 
group sense, i.e.: further contributions to as well as non-Gaussian and non-Markovian 
parts of the Langevin forces are irrelevant for the long-time and large-distance properties of 
the system as long as the diffusion constants are positive. 

For an infinite system the term proportional to the coupling constant / in the Langevin 
equation (^) can be eliminated by a suitable Galilean transformation ll^]. Such a Galilean 



transformation can of course not be applied to a system with time independent surfaces 
perpendicular to the driving force. However, Krug has shown |jl3| that in a system with 
open boundaries in which the particles are driven from a reservoir of density ci to a second 
reservoir of density C2 with ci > c* > C2 the particle density far in the bulk takes the 
value c* which maximizes the function k(c). In this case we have / oc — /t'(c*) = and 
g oc — k"(c*) > 0, and the density decreases only algebraically from the boundary value 
to the bulk value. Since we are interested in the behavior of the system near one of the 
boundaries we may effectively consider a semi-infinite system with bulk density c = c*, i.e. 

/ = o. 

Appropriate boundary conditions for diffusive semi-infinite systems have been derived 
in [0,1). This approach crucially rests on the assumption of detailed balance. In Ref. |Tl| 
it was shown that the dynamics defined by the bulk stochastic equations of motion 
with / = satisfies this assumption ifw = cr/p = l. Under renormalization group transfor- 
mations the coupling constant w flows to the detailed balance value = 1 ii d < 2. For 
d > 2 the system displays mean field-like behavior. It should be remarked that for w 1 
the correlations are long-ranged ||19|. Thus for d < 2 the variable {w — 1) plays the role of 
a dangerously irrelevant coupling and its effects have to be studied. 

To set up a renormalized field theory it is convenient to recast the bulk model in terms 



of the dynamic functional pO|-p5|JTT 



Jb[s,s\ = J dtJ^d\[sdtS + \[iV±s)(y±s) 

+p(dfs)idiis) + fidfs)s + ^gidfs)s' 



{Vrs)''-cx{dfs)' 



(7) 



where s is a Martin-Siggia-Rose response field [^. The range of integration with respect 
to r is the ci-dimensional half-space ^ = {r = {r±,z) \ r± G R"'^^,^ > 0}. Correlation 
and response functions in the bulk can now be expressed as functional averages with weight 
exp(— jTi,). For w = a/ p = 1 the functional jTj, can be written in the detailed balance 
form [g4|,g5 

6H 



dt d'^r [s{dts + - s 

with the reaction kernel 

R, = A[(V±V± + p\d\\) + |(9„s - sd\\)]. 



(8) 



(9) 



Here V and d act to the left, while V and d act as usual to the right. The Hamiltonian 

H[s] = j^d'^r\s{v)\ (10) 

defines the purely Gaussian stationary state distribution. Obviously, jTj, now obeys the 
detailed balance symmetry P^j2^ 



s{t) -s{-t) 

m s{-t) - ^ 

OS 



si-t)-si-t). 



(11) 



The last equation implies in particular that 



{siv,t)§iv',t')) = Bit -f){siv,t)siv',t')). 



(12) 



We now turn to the surface contributions J'l in the full dynamic functional J = Jh^ J\ 
arising from the boundary layer at 2; = 0. Locality is assumed for the bulk and the surface. 



Therefore J\ should be written as a (ci — l)-dimensional integral over surface fields alone. It 
is easily seen that for the Hamiltonian 7i only irrelevant surface terms can be constructed. 
Thus 7i has the form (^) also in the semi-infinite case. To retain the detailed balance 
form (^) of J we can only modify the reaction kernel (y) R;, ^ R = R;, + Ri. Such a 
modification has to describe the breakdown of particle conservation due to the reservoir . 
The relevant contributions which respect the symmetry (|ll|) are 

Ri = A(c + 6(s -25))5(z) (13) 

up to redundant terms. Since the kernel R has to be positive for reasons of stability it 
follows that c > 0. This completes our construction of the dynamic functional in the case 
of detailed balance. We eventually get 



I dt l^d''r[sdts + x{iV±s){V±s) + p(9||5)(a||s)) + ^Xg{diis)s^ - A((Vx5)' + pidfs, 
+ f dt f^^d'^~^r±x[cs{s - s) + bs{s - s){s -2s) -^gss^]. (14) 



dV 

We now consider the modifications that we expect if detailed balance does not hold. First 
of all, the bulk functional is given by (^ with a noise constant a independent of the dif- 
fusion constant p. In addition, the different surface fields which show up in ( p^ ) are now 
independent, and a boundary source has to be added to the equation of motion (|) leading 
to a term linear in si. We thus get for the surface dynamic functional 



J7i[s, s] = J dt J d'^ ^r±X(^css — cs^ — his 



+ ]j^9aSS^ - ]^gbS^s + ^9cS^) ■ (15) 



In Eqs. ([T^,|T^) we have omitted the surface operators S(9„s, (9„S, and d'^s (where 
dn = <9|| means the normal derivative) since they can be expressed in terms of the composite 
fields retained in J7i (see also below). 
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We remark that we have always used the prepoint time discretization in the construction 
of the dynamic functional p^ , p5| . This corresponds to the definition Q{t = 0) = and thus 
allows to omit all (measure) terms oc 0(0) in J'. 



III. EQUATIONS OF MOTION AND MEAN FIELD PROFILE 

Introducing the functional 

Z[J, J;Ji,Ji] 

= J Vis, s] exp(^-Jb[s, s] - Ji[ss, Ss] 

+ ( J, s) + (J, s) + ( Ji, Ss) + ( Jl, Ss)) (16) 

we may write correlation and response functions as derivatives of Z with respect to the bulk 
sources J, J and the surface sources Ji, Ji. In Eq. (|1^), Sg and Sg denote the surface fields 
and we have used the abbreviations 

(J,s) = J dtJ^d'^rJS (17) 

and 

{J,,Ss) = [dt [ d^-'r^J,Ss. (18) 
J JdV 

To obtain the boundary conditions which are determined by the surface functional J^i we 
exploit the invariance of Z with respect to a shift of the fields s and s. This invariance 
implies the equations of motion 

{^Ja){jj.j^j^} = Ja for A = s, s, Ss, Ss, (19) 

with the notation Jg = J , Jg = J , Jss = Ji and Jg^ = Ji respectively, and 

6Js = dts - A((Ax + pd^)s + fd\is + ^gd^ys^ 

-2{A^ + adf^s), (20a) 
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5 = -dtS - A (A^ + pdXs - fdfs - g{d\{s)s), (20b) 



= \ {^pdnS + 2adnS + {c- f)ss - 2css 
+^{9a - g)sl - gbSsSs + ^gJl - hi) , (20c) 

SJs, = A (-pi9„s + cSs + gaSsSs - ^gbsf) ■ (20d) 



Equations ([l9| , |20c| . |20dD show that dnS and dnS can be expressed in terms of the surface 



operators included in J'h. The redundance of SgdnS, SgdnS, and SsdnS follows by differentiating 
{dns) and {dns) with respect to Ji and Ji, respectively. 

Setting J = J = Ji = Ji = in ([T9|) we obtain from (|20a|) the equation 



d,{pdM^) + m^) + l9{s{r±,z;tf))=0 (21) 

for the steady density profile ^{z) = {s{r±, z;t)) . Using the boundary condition $buik = 
which follows from our definition of s an integration with respect to z yields 

p^'iz) + mz) + lg{Ciz) - Cbuik + Hzf) = 0, (22) 

where the function C{z) = ((s(r_L, z; t) — $(z))^) (with the bulk value Cbuik) describes density 
fluctuations. The boundary condition at 2; = follows from ( p^Ocj ) and (^) as 



C$0 + lga{CiO) + - i^?Cbuik = hi. (23) 

Note that equal time expectation values such as (s(t)"5(t)™') with m > are zero as a 
consequence of our prepoint time discretization. The function C{z) contains an ultraviolet 
divergence proportional to 5(0) = A^^^Ay (where A_l and Ay are cut-off wave numbers) 
which cancels in Eq. (^2]) for z > 0. 

The mean field profile is the solution of Eq. ( P^ for vanishing fluctuations, i.e. C{z) = 0. 
In the case / = which describes the maximum current phase we obtain 

$^f(z) = $o(l + |-$o^)"'. (24) 
zp 
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IV. RENORMALIZATION GROUP ANALYSIS 



A. Surface divergencies 

In the next subsection the density profile will be calculated by an expansion around 
the mean field profile (0). Since the perturbative calculation of Green functions leads to 
integrals which are ultraviolet-divergent in the upper critical dimension = 2, we have to 
employ a regularization method to obtain a well-defined perturbation series. We utilize the 
method of dimensional regularization which allows us to render all integrals finite up to poles 
in e = 2 — d. The e-poles are then absorbed into renormalizations of coupling coefficients 



and fields (minimal renormalization) . It was shown in Ref. |TT| that for the bulk system a 
renormalization of the coupling coefficients p and a is sufficient to cancel all uv-divergencies. 

In general the breaking of translational invariance by a boundary generates additional 
divergencies which require additional renormalizations of surface couplings and fields. (For 
a review see Ref. [|I|.) To determine the surface divergencies we first consider the fiat profile 
^{z) = which is a steady state of the system provided that $o = and C{z) = Cbuik for 
all 2; > [see Eq. ([2^)]. [At the stable fixed point with w = a/ p = 1 the z-independence of 
C{z) follows immediately from the the Hamiltonian (|T^).] It will be shown below that this 
condition is generally satisfied for d = 1 and in the high-temperature, strong-field limit in 
higher dimensions. 

The Fourier transform G'q^^^(z, z') of the free (Gaussian) propagator (s(r_L, z; t)s{0, z'; 0)) 
is the solution of the differential equation 

{tcu + Ml - Ap92)Gq^,^(^, z') = 5{z - z') (25) 

with the boundary condition 

pdz'Gci^,uj{z, z) ^^^^ = cG'q^,^^(z, 0). (26) 

This follows in a straightforward way from the terms bilinear in s and s of the functional 



(with / = 0) and Eqs. ( p!^ , pOd|) . The solution is given by 
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Gc,^,^{z,z') (27) 



2\^K L K + c/ ^ 

with 

/iClJ r,\l/2 ^ ^ 

« = (y + ?i) • (28) 

The parameter c occurring in the surface functional Ji and in the propagator describes (for 
c > 0) the suppression of density fluctuations by the particle reservoir at the boundary 
z = Q. Since its momentum dimension is c ~ ^/p^^ the asymptotic scaling behavior is 
governed by the fixed point = oo. This means that up to corrections to scaling the 
boundary condition (p6|) may be replaced by the Dirichlet boundary conditions s{z = 0) = 
and s{z = 0) = 0. 

The Fourier transform of the Gaussian correlator (s(r_L, t)s(0, 2;'; 0)) at the Dirichlet 
fixed point follows from (|^) as 

Cc,^A^,z') (29) 
= 2A / dyGq^^^{z, y){q]_ + adydy)G-^^-^{z\ y). 

J 



Using the differential equation (|25|) for the propagator and the Dirichlet boundary conditions 
this may be written in the form 

C'q^,.(^, Z') = W[G^^,^{Z, Z') + G'-q^,-.(^', z)] (30) 

+2(1 - w)\q\ / dyGq^,^{z, y)G-^^,-^{z, y), 

•J 

where w = a/ p. For the calculation of density fluctuations it is useful to transform the 
correlator into (q_L, )f:)-representation. The result reads 

C^A^,z';t) (31) 

= wG^Az,z'; \t\) + {l-w)\ql / G^^z, z' ;t')dt' , 

1^1 



where 
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xe 



iz-z')\ / (z + z')^ 



exp I — I — exp 



AXpt ' V 4Apt 

is the free progagator. 

For w = \ Eq. (|3TD reflects the detailed balance symmetry (|12|). Integrating the correla- 
tor (|3T|) over we obtain for t = 0, z > 

- = -(1 - ^''^ 
It is thus necessary for the flat profile to be stationary in dimensions d> \ that the condition 
= 1 is satisfied. Deviations from this stable fixed point generate the long-range correla- 
tions IP] in the bulk which are responsible for the z-dependence of the function C(z). In this 
sense (w — 1) is dangerously irrelevant. There is, however, a microscopic realisation of DDS 
in which equal-time correlations are absent in spite of the violation of (microscopic) detailed 
balance p^: In the high-temperature, strong- field limit the particles are non-interacting ex- 



cept for the excluded volume constraint and cannot jump against the direction of the driving 
field. The lack of correlations in this limit can be incorporated into the field-theoretic de- 
scription by setting w = \ from the outset. For ci = 1 we can render w = 1 by a simple 
rescaling of coupling coefficients and fields. 

We now proceed along the lines taken in the analysis of the surface critical behavior of 
equilibrium systems at the ordinary transition P|J^,P^. Since the correlator ( pi] ) and the 



propagator ( P^ vanish for 2;' = or 2; = 0, Green functions with insertions of the surface 
fields Sg and Sg tend to zero for c ^ 00. A convenient method to investigate the scaling 
behavior of quantities that vanish for c = 00 is the 1 / c-expansion . 
For large c the propagator of the surface field behaves as 

G^,^^{z,^)=c-^ pd,G^^lS^,A\,^,^0{c-^) (34) 

for z > 0. In order to derive an analogous relation for the correlator we have to take into 
account the surface coupling c in J\ which has the same //-dimension as c. By differentiating 



the equation of motion (|^, |20c| ) with respect to J one obtains 
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z'=0 



2a d^>Gci^^^{z,z'] 



z'=0 



(35) 



+ cC^^^^{z,0) -2dG^^^^{z,0) 



and, for c — >• cxD, 



C^,,M,0) = c~'p[d.,C^^%{z,z') 



z'=0 



-2iw-w,)d,,Gi^2^{z,z' 



z'=0. 



+ 0(c- 



(36) 



(with Ws = c/c). Therefore the leading order terms in an expansion in powers of c"^ can 
be studied in the framework of a field theory with Dirichlet boundary conditions by the 
replacements 



C ^pdnS 

Ss c~^p{dnS - 2(-u; - Ws)dnS) 



(37a) 
(37b) 



in expectation values. To investigate the scaling behavior of Green functions with insertions 
of Ss and Ss we have to compute the uv-singularities generated by the surface fields on the 
r.h.s. of ([37|). These singularities require renormalizations of the form 



pdnS = Z\''^[pdnS\R 

p{dnS - 2wdns) = Z\''^[p{dns - 2w(9„s)]r 

^1/2 ryl/2 



(38a) 
(38b) 
(38c) 



The renormalization factors Zi and Zi are determined by requiring that the averages 



{[pdn5]R-s{z > 0)) 



and 



{[pidnS - 2wdns)]R ■ s{z > 0)) 

be finite for e = 0. Due to causality the term 2pwdnS gives no contribution to the second 
expectation value but in the Green function {p{dnS — 2wdns) ■ s{z > 0)) it is necessary to 
obtain a surface field that can be multiplicatively renormalized by Zi. 
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The renormalizations which are necessary to cancel all uv- divergencies in the transla- 
tionally invariant bulk theory are given by 



p —>■ p = Zpp a —>■ a = Z^o (39) 



with 



u 



= 1 - + + 0(^2) (40a) 
8e 
u 



where 



^^■(3(w; + w;-1) + 2)+0(m2) (40b) 
loe 



In equation ( [40c| ) , p is an external momentum scale and = r(l — e/2)/((l + e)(47r)'^/^) is a 
geometrical factor which has been introduced for convenience. The renormalization factors 
Zi and Zi are calculated at one-loop order in Appendix with result 

^1 = 1 + ^^^ + 0(^2) Zi = 1 + 0(^2). (41) 

We now show that at the stable fixed point w = 1 the theory is free of surface singularities, 
i.e. Zi = Zi = 1. The idea for the proof is presented graphically in Fig. |l|. The hatched 
bubble represents the sum of all connected Feynman graphs with a single external s-leg and 
two s-legs. Together with the three-point vertex representing the bulk coupling {Xg/2){d\\s)s'^ 
it gives the sum of all perturbative contributions to the (amputated) Green function Gi_i. 
Lines with or without an arrow denote Dirichlet propagators or correlators while a short line 
perpendicular to a propagator /correlator line indicates a derivative with respect to z. By 
naive power counting one expects that the integrations associated with the first Feynman 
graph diverge as Ay. Since the breaking of translational invariance reduces the degree 
of divergency by unity ||1[] the graph should give rise to logarithmic surface divergencies. 
However, Fig. |I] shows that the graph may be replaced by a different diagram with only 
logarithmic bulk singularities (and thus without surface divergencies). 
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The first equality in Fig. ^ is justified if we assume that the time variable carried by the 
leftmost vertex is smaller than the time variables associated with the vertices in the bubble. 
In this case we can replace the first correlator by a propagator [Eq. (|1^)]. In fact, if we 
assume that an internal vertex carries the earliest time argument the graph vanishes due to 
the Dirichlet boundary conditions (see Fig. ^). The second equality in Fig. |^ follows from 
integration by parts with respect to the 2;-variable of the first vertex: 

/ dzssdzS = / dzs-dz{s^) = — - / dzs'^dz-s, 
Jo Jo 2 2 JO 

where we have used the Dirichlet boundary conditions. Since the 2;-derivative now acts on 
an external line the degree of divergency is reduced by unity. The remaining logarithmic 
bulk divergency is cancelled by the renormalization of p. 

B. Density profile at one-loop order 

In order to compute the lowest order correction to the mean field profile we insert the 
ansatz 

^z) = ^^i{z) + ^l{z) (42) 
into Eq. (^) and retain only the terms linear in ^i{z). The resulting differential equation 

p^[{z) + ^9{C{z) - Cbuik + 2^^i{z)^i{z)) = (43) 

has the solution 

$i(z) 

= -TT dz'expi — / ^miiy)dy]{C{z') - Chnik) 
2p Jo ^ p Jz' ' 

where a = 2p/{g^o). At one-loop order it is sufficient to calculate the density fluctuation 
C{z) in a Gaussian approximation. Since this function depends on the density profile we 
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have to generalize the results of the previous subsection to the case $(-2) 7^ 0. For this 
purpose we split the field s(r) into its mean value and a fiuctuating part 0(r), which 
leads to the dynamic functional 



(45) 



where 



(46) 



cit / 

V 



<pdt<P + A((9x0)(9x0) + p{d\\<P){d\\<P) 
+g^[z){d{^)<P-{d^~^f-p{d{^f), 



(47) 



and 



int ) ' 



(48) 



(Throughout this subsection we assume w = 1, i.e. a = p.) The Gaussian propagator 
satisfies the differential equations 



[k' - pdl - gdMz)]G^^A^. z') = -5{z - z') 
[k" - pdl + g^{z)d,]G^,,^{z' , z) = ^5{z - z') 



(49a) 
(49b) 

which follow from the functional Jg (^TI). The solution for $(-2) = $mf(-2) may be written 
in the form 



1 1 



(50) 



/+(Co 
/-(Co 



-/4C<)-/+(C<) 



where 



/±(C) = (lTC)exp(±C) 



(51) 
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and 



C = K{z + a)/^, (Q = Ka/^, (52a) 
= K(max{z, z'} + a) / y/p, (52b) 
C< = K{mm{z, z'} + a)/ ^/p. (52c) 



To obtain the correlator one has to calculate the integral 



Cci^,co{z,z') =2A / dyG^^^^{z,y) 

•J 

x{ql + pdydy)G^^^.^{z\y). (53) 



We use Eqs. (|^) to simplify the integral by eliminating one of the two ^/-derivatives in 
The result 

C^^,^{z, z') = G^^^^^iz', z) + G^^^^{z, z') (54) 

POO 

-Xg / dy^^i{y)dy[Gq^^^{z,y)G^^^_^{z',y)] 

J 

can be further simplified by integration by parts since 

G^^,_^{z\y)dy^^,{y) = Gq„-a.(l/,^')'9.'$mf(^')- (55) 
This manipulation leads to 

Co,^,u,{z, z) = Gq^ -c^(2;', z) + Gci^^^{z, z) (56) 
2Ap 



' / I ^2 / <i'yGc,^,u,{z,y)G^^^-^[y,z!). 

The function C(z) describing the Gaussian density fluctuations is given by the integral of 
C'q^,a;(-2, -z) ovcr q_L and ijj. Using the semi-group property 



Gqx(^) Z ] t\ + t^) 

= / dyGq^{z,y;ti)G^^{y, z';t2) (57) 
of the propagator in (q_L, 2;, t) -represent at ion we obtain 
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C{z) — Cbuik — 



\p 



z + a) 



(58) 



where 



(59) 



is independent of 2; > 0. The straightforward evaluation of the g^-integral in (^) yields 

„ r(i-6/2)2p(^-)/^ ,r(^ + «)' 

o — Obuik = . ^ ./.. — ^ —z 



az 



(47r)'^/2 (;z + ay 



dy- 



l + e/2 r(l-e)7o ^l + ay/(22 
To compute the one-loop correction to the profile one has to insert (60) into Eq. 

r(l-e/2) 9l2-y-'/' ^zjfxz/^y 



(471)^^/2 (z + ay V e(l + e) 



— F{z; a, e 



with 



F{z; a, e) = (/i/ ^y / dz' 

'0 (z +a)^ 



/ 1 2 

(— + 



+ e/2 r(l-e)yo y + 2z'/a 



y - 

dy _ . , e ^ 



In the limit e ^ we have 



(60) 
giving 
(61) 



(62) 



F{z;a,e = 0) = z(l + 



z + a 

/ , ,22;, z — a ,2z.^ ,2z. 
-ah + In — - —— exp — El — 
^ a z + a a a 



(63) 



where 7 = 0.5772... is Euler's constant and Ei(x) denotes the exponential integral 



Eq. ( [60| ) shows that the perturbation series is not well-defined for e = 0. This uv divergency 
is cancelled by the renormalization of the diffusion constant p in the mean field profile which 
leads to 



-iif^z/Vpy-i] 



2p 1 r u 
g z + a^ 2[z + aj ^ e 
-(1 + e)F{z; a, e))] + 0(2 - loop) 



(64) 
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Since the coupling constant g is relevant below two dimensions we may not use the per- 
turbation series at finite order directly to study the asymptotic behavior of the profile. In 
the following subsection we improve the one-loop result Eq (Q) by a renormalization group 
analysis. 



C. Scaling 

The renormalization group transformation allows us to map the profile for large values 
of z (compared to microscopic length scales) to length scales which are accessible to pertur- 
bative calculations. Since the bare coupling coefficients are independent of the momentum 
scale /i the profile satisfies (for w = 1) the renormalization group equation 



[nd^ + P{u)du + Ciu){pdp + ada)]^{z;u, p, a; jj) = 0, 



(65) 



where ||Tl| 



C{u) = - p-r- 
p dp 



P = --u + 0{u' 



and 



p{u) = p 



d_ 

dp 



u = u{-e--C{u)) 



(66) 



(67) 



are defined as derivatives at fixed bare parameters. Eq. (|65|) can be solved by the method 
of characteristics, with the solution 



$(z; M, p, a; p) = ^{z; u{l) , pX (l) , aX{l); pi) 



(68) 



and associated characteristics 



(69) 



Using dimensional analysis and Eq. (^) it is easy to show that the profile satisfies the 
relation 
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^{z■u,p,a■^^) = {^llYI\pX{l))-'/' (70) 
x$(/i/z(pX(/))-i/2. iila{X{l)/pY'^- 1). 

For small / the scale dependent coupling coefficient u{l) tends to the fixed point Ui, = 
(4/3)e + O(e^) and the function X{1) is asymptotically proportional to a power of /, namely 

x{i) ^ Xor'\ V = -^C(^.) = (71) 

where Xq is a non-universal scale factor. By choosing the value 

/ = [iXopY/Vipz)f^'-'^ (72) 
for the fiow parameter we obtain for z ^ \fpp~^ the scaling form 

$0) = $oH(5<l>oz(''+')/('-'')) (73) 

with the scaling function 

H(x) = x"^$(1;m^, l,aa;"^; 1) (74) 
(where a = [^A^jUi)^!'^') and the non-universal constant 

i? = ((pXo)i/V/')''^'"'^^^'^'"'"- (75) 
To derive Eq. (|7^) we have used the identity 

A^g" = u,{pXof/^p^ (76) 

which follows from the fiow equations (|69D and the initial condition ( |40c|) to express a in 
terms of $05 Pi 1^*, Xq, and p. Note that for d = 1 there is no adjustable parameter in 
Eq. (|73| ) since B is unity in this case. The scaling behavior of the one-dimensional profile 
was first predicted by Krug |T3|| . 

The reciprocal of the scaling function at first order in e is given by 

H(x)-i = 1 + - (1 + ^F(l; a/x, e = 0)) + ©(e^). (77) 
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In order to compare this result with the exact solution found in Refs. [14,15| we have to take 



the continuum limit of the exact profile. In Appendix ^ we show that this limit leads to 

2exact(a;) = exp(a;^)erfc(a;), (78) 

where erfc(x) = J!^ dy expl—y"^) denotes the error function [^. The exact scaling func- 
tion and the one-loop approximation ([77|) for e = 1 are depicted in Fig. ^. The asymptotic 
form 



■^exact 



X] 



^ ^/nx hi X ^ oo (79) 



of the scaling function can be compared with the asymptotic behavior 

E{x)-^ ~ -(1 + + 0{e'^))x ^ l.O860Fa; (80) 
a 3 

of the e-expansion at first order for e = 1. Here and in Fig. ^ we have used the geometric 
prefactor a = (474^/^^,)^/^ directly for e = 1 (i.e., without expansion in e) since this procedure 
is empirically known to give the best numerical results [pT|. 

By a different choice of the flow parameter it is possible to eliminate the 0(e)-contribution 
to the profile completely. The corresponding value of / is defined by the equation [see ( p^ , 

m] 



IJ,lz{pX{l))-^/^ \n{filz{pX{l))-^/^) (81) 
-F[filz{pX{l))~^/^; pla{X{l)/ pY^^ ,e = o) = 0. 

However, this modification has only a small effect on the scaling function ^(s)"^ at one-loop 
order. 

We conclude this section with the discussion of logarithmic corrections to the mean field 
profile in two dimensions. For e = 0, w = 1 the scale dependent coupling coefficient u{l) is 
the solution of the differential equation 
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with the initial condition u{l) = u. For / — > we find the asymptotic expression 



u 



(0^i(.n(l/0r'[l+O(lil|i^)], (83) 



A straightforward calculation yields the asymptotic behavior of the characteristic X{1): 

lnX(0 = / -yC,{u{l')) ~ -lnln(l//) + cst (84) 

for / 0, i.e. X{1) ~ (ln(l//))2/3. In order to sum the logarithmic singularities in the 
perturbation series we set I = ^/p/ (fiz) in Eq. (|68| ) and obtain 

^z, $o) = $0 [l + est X $oz(ln(/i^/^))-'/'] ' (85) 

(lnz)2/3 ^ ^ ^ 

~ ^ — for z^oo. (86) 

This result is valid on length scales that are large compared to microscopic distances, i.e. 

z^A^K 

V. SUMMARY AND OUTLOOK 

Until now most analytical studies of driven diffusion deal with closed systems. To main- 
tain a steady state with a non-vanishing current one usually imposes periodic boundary 
conditions in the direction parallel to the driving force |jlO|- Since in real physical systems 



the particles are fed into the system by a reservoir the effects of boundaries are of great 
interest. Using a field theoretic renormalization group approach we have calculated the 
density profile in a semi-infinite DDS with a particle reservoir at the boundary. For d = 1 
our result is in fair agreement with the exact solution for the totally asymmetric exclusion 
model [l^JT^. Equation (|85|) for the two-dimensional profile is a new result which may be 



checked by computer simulations. 

We plan to extend the present work along several lines: 

1) We have studied the effect of a particle source (or sink) at z = assuming a second 
reservoir to be located at z = L with L — oo. One would also like to consider the interplay 



of the reservoirs for finite L \14A5 
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2) In a system with quenched disorder and periodic boundary conditions a particle is 
subjected to the same random potential after every passage through the system. In this 
way the periodicity produces long-range correlations which make it difficult to compare 



simulational results with theoretical predictions which assume uncorrelated disorder |]T7 . 
To circumvent this problem it would be desirable to consider disordered DDS with open 
boundary conditions. 

3) In the case of non-critical DDS the scaling behavior (|73|) of the proffie can be described 
by a single exponent, namely the bulk exponent rj = {2—d)/3. This is due to the fact that the 
surface density $o has the same scaling dimension as the bulk density. We expect that this 
is no longer true if the system is at its critical point. In general the surface critical behavior 
is governed by new exponents which cannot be expressed in terms of bulk exponents. 

4) Another line of extension would be to consider boundaries parallel to the driving force. 



5) In Sect. |^ we have shown that for d > 1 the flat profile ^{z) =0 is only stationary 
in the case of detailed balance and in the limit of high temperature and strong field. In 
the more general case of a lattice gas at finite temperature the longe-range correlations in 
the bulk give rise to an inhomogeneous steady state even for $o = $buik- (It was already 
pointed out by Zia and Schmittmann |]TU[ that the coupling coefficient {1 — w) which is 
responsible for the long-range correlations is dangerously irrelevant for 1 < c/ < 2.) This 
effect should lead to a modification of the maximum current principle [Eq. (|lD]. Finally, it 
would be interesting to study the formation of the stationary density profile starting from a 
homogeneous initial state. 
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APPENDIX A: SURFACE SINGULARITIES AT ONE-LOOP ORDER 



In Sect. |T^ we have shown that for w = 1 the perturbation expansion is free of surface 
singularities. In this appendix we consider the general case w > 0. 
The Fourier transform of the surface-bulk response function 

Xsb{r±, z; t) = (s(r_L, z; t) ■ p(9„5(0, 0)) (Al) 
and the bulk-surface resonse function 

Xbsir±, z; t) = {pdnsir^, t) ■ 5(0, z; 0)) (A2) 
may be written in the form 



z'=Q 



Xsb{q±, u;,z) = p d^'G^^^^{z, z) 

roo poo 

+ dy dy'K^^,^{y,y')dyGq^,^{z,y) pd,'Gq^,^{y',z') , (A3) 

JQ Jo ^ — ^ 



and 



z'=0 



X6s(q±, c^, z) = p d^>Gq^^^{z', z) 
+ dy dy'K^^,^{y',y)dyG^^,^{y,z) pd^>G^^^^{z',y') (A4) 

JO Jo 2 =0 

respectively, where Gq^^i^{z, z') is the free propagator (p2| ) in Fourier representation. At 
one-loop order the kernel Kq^^^{z, z') is given by 



') = {\gf I [9,,G,,_,,,^_.(z, z')]C'k„.(^, z) + 0(/). (A5) 



Note that K^^^^^{z, z') has to be treated as a distribution since the calculation of response 
functions involves integrations of Kc^^^^^{z, z') with respect to z and z' . 

Before we calculate the singular part of the kernel explicitly it is helpful to discuss the 
general form of the divergencies one has to expect. First, K^^^i^{z, z') has a non-integrable 
(for d = 2) singularity if z —>■ z'. This divergency occurs in the translationally invariant bulk 
theory as well as in the semi-infinite model and has to be cured by a renormalization of p. 
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The form of this singularity follows from dimensional analysis and the isotropy in the direc- 
tions parallel to the surface. Since K(^^^^{z, z') ~ AyU^ and z ~ ^/pfJ'^^ the (dimensionally 
regularized) bulk singularity contained in K^^^^^z, z') is proportional to \p5\z — z'). The 
breaking of translational invariance by the boundary generates additional divergencies for 
z, z' ^ {) [|1|]. Upon dimensional regularization they are proportional to \p5{z)5{z'). 

In order to compute the coefficients of the distributions we apply the kernel to exponential 
test functions. At lowest order in g"^ we have 



dz dz'K„^^{z,z')e-^'-^' 
Jo 



V Ac/Vr(l + e/2) p3 + w (3' 



(A6) 



The inverse Laplace transform of this expression with respect to (3 and (3' is given by 



A^Vr(l + e/2) r3 + w.,, , w-1 



where we have introduced the definition 



S'{z'\z) + ^6{z)6{z') +0{e',g'), (A7) 



dz'S'{z'\z)f{z') = ~f\z). 



(A8) 



(This distribution is equivalent to 5'{z' — z) if we define S{—z) = for z > 0.) It is now 
straightforward to calculate the singular parts of the response functions at one-loop order. 
Inserting (^) into Eqs. ([A3|) and (|A^) we find 



Xsb{ci±,uj,z) = -^1 + -(^ 



1 r U /-W — 1 3 + W KZ 

^7*^ 8 ^ 16 7^ 



--/VP 



(A9) 



and 



Xbs[fl±,uJ,z) = - 



(AlO) 



e 16 

To obtain renormalized response functions one has to express the coupling coefficients p and 
a by their renormalized counterparts [see Eq. (|39|)] and absorb the remaining e-poles into 
renormalizations of surface fields, i.e. 



[Xsb]R — Zi ^^"^Xsh, 



[Xbs]R — ^^'^Xbs 



(All) 
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with 

~ 111 — 1 ?/ 

Zi = l + ^^- + 0{u^), Zi = 1 + 0(^2). (A12) 

The renormahzation factor Zi is unity at every order of the perturbation series since the 
primitive surface divergencies in Kf^^^i^{z, z') are proportional to 6{z)6{z'). Due to the Dirich- 
let boundary conditions this distribution has no effect in Eq. 



APPENDIX B: THE EXACT PROFILE IN THE CONTINUUM LIMIT 

One of the simplest examples of a DDS is the totally asymmetric exclusion model in one 
dimension . For this model with nearest neighbor hopping and open boundary conditions 
the density profile has been calculated exactly by Schiitz and Domany |]T3[ and Derrida et 

ai. m 



The model is defined on a chain of sites each of which is either occupied by a single 
particle or empty. At each time step one chooses a random lattice site p, 1 < p < N. If this 
site is occupied by a particle {up = 1) and its right neighbor site is empty (rip+i = 0), the 
particle will jump from p to p + 1. In the model with open boundary conditions particles 
are injected with rate a at the left boundary and leave the system with rate /5 at the right 
boundary. 

For > 1/2 the bulk density takes the critical value ric = 1/2 and the system carries 
the maximum current. Since we wish to compare the e-expansion of Sect. |V|with the exact 
profile we let A^ oo in order to eliminate the infiuence of the right boundary. Using the 
formulas given in the appendix of Ref. we arrive at the particle density 



1 ^,3 ^ 4a(l-a) 



2p(2a-l)2 ^ '2'^ ' (2«-l)2 
for a, /5 > 1/2 and A^ oo, where 6; c; z) is a hypergeometric function 

To compute the scaling function Sexact(a;) from the microscopic profile one has to take an 
appropriate continuum limit. In addition to the distance from the boundary the injection 
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rate a provides the only macroscopic length scale = (2q; — in the model. We 
therefore measure p in units of and let p, — » oo at fixed = p/^a- For this purpose it 
is convenient to use the integral representation 

=p f\i-tr-\i + {^^-i)t)-'/'dt (B2) 

Jo 
fp 

= / exp((p — 1) ln(l — r/p)) 
Jo 

x(l + (e„-l)r/p)-3/'rfr (B3) 



to show that 



oo 

e-"(l + r/x2)-3/2rfr. (B4) 



For large p the prefactor in Eq. ( pTD is given by 



A straightforward calculation now yields in the continuum limit defined above 

(o-p) = (Toexp((T^p)erfc((Jov^), (B6) 

where we have expressed the occupation numbers in terms of the spin variables ap = 2np — 1, 
(To = 2a; — 1. 
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FIGURES 





FIG. 2. The vertex is assumed to earry the earhest time argument. In this case the integration 
over z vanishes due to the Dirichlet boundary conditions. 
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(a) (b) 



FIG. 4. Perturbative contributions to the (a) surface-bulk and the (b) bulk-surface response 
function. Here the hatched bubble represents the integral kernel (|A5|). 
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